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ABSTRACT 



lO \ We investigate the effects of a passing stellar encounter on a planetesimal disk 

through analytical calculations and numerical simulations, and derive the boundary 
^ ' radius (opianet) outside which planet formation is inhibited by disruptive collisions with 

high relative velocities. Ida, Larwood, and Burkert (2000. ApJ. 528, 1013-1025) sug- 



^^ ' gested that a stellar encounter caused inhibition of planet formation in the outer part of 



a protoplanetary disk. We study orbital eccentricity (e) and inclination (i) of planetes- 
imals pumped up by perturbations of a passing single star. We also study the degree of 
alignment of longitude of pericenter and ascending node to estimate relative velocities 
between the planetesimals. We model a protoplanetary system as a disk of massless 
^ particles circularly orbiting a host star, following Ida et al. (2000). The massless par- 

\-{ ' tides represent planetesimals. A single star as massive as the host star encounters the 

^ ■ protoplanetary system. Numerical orbital simulations show that in the inner region 

at semimajor axis a < Q.2D where D is pericenter distance of the encounter, e and 
i have power-law dependence on (a/D) as e oc (a/D)^''^ and i oc (a/D)^''^ and the 
longitudes are aligned, independent of the encounter parameters. In the outer region 
5^ I a > 0.2D, the radial gradient is steeper, and is not expressed by a single power-law. 

The longitudes are not aligned. Since planet accretion is inhibited by e as small as 0.01, 
we focus on the weakly perturbed inner region. We analytically reproduce the power- 
law dependence and explicitly give numerical factors of the power-law dependence as 
functions of encounter parameters. We derive the boundary radius (cpianet) of planet 
forming region as a function of dynamical parameters of a stellar cluster, assuming the 
protoplanetary system belongs to the stellar cluster. Since the radial gradient of e is 
so steep that the boundary is sharply determined. Planetesimal orbits are significantly 
modified beyond the boundary, while they are almost intact inside the boundary. This 
tendency is strengthened by reduction of relative velocity due to the longitude align- 
ment in the inner region. We find apianet ~ 40-60AU in the case of Z? ~ 150-200AU. 
D ~ 200AU may be likely to occur in a relatively dense cluster. We point out that the 
size of planetary systems (flpianet) born in a dense cluster may be necessarily restricted 
to that comparable to the size of planet region (~ 30-40AU) of our Solar system. 



Subject headings: celestial mechanics; orbits; planetesinial; accretion 



1. INTRODUCTION 

In general, stars are born as members of an open cluster. Stellar clusters would evaporate 
on timescales more than 10® years (Kroupa 1995; 1998). This evaporation would be caused by 
gravitational interactions between stars, so that many stars experience gravitational perturbations 
of the other stars during the evaporation. More than half of T Tauri stars have protoplanetary disks 
(e.g., Beckwith and Sargent 1996), which would eventually form planetary systems on timescales 
10^-10^ years (e.g., Safronov 1969; Wetherill 1980; Hayashi et al. 1985). Planetary systems would 
be affected by stellar encounters more or less during their formation stage. 

In the standard model (e.g., Safronov 1969; Wetherill 1980; Hayashi et al. 1985; Lissauer and 
Stewart 1993), terrestrial planets and cores of jovian planets accrete from planetesimals that are 
formed in a protoplanetary disk. The accretion of jovian planet cores is followed by gas accretion 
onto the core when the core acquires critical mass ~ 5-15 M^ (e.g., Mizuno 1980; Bodenheimer 
and Pollack 1986; Ikoma et al. 2000). 

The passing stellar encounters would pump up orbital eccentricity e and inclination i of plan- 
etesimals. The velocity dispersion of planetesimals is given by ~ Ve^ -|- i^fkep; where i is given in 
unit of radian and i;kep is Keplerian velocity (e.g., Safronov 1969; Lissauer and Stewart 1993; Oht- 
suki et al. 1993). If the velocity dispersion exceeds their surface escape velocity of planetesimals, a 
collision between the planetesimals results in disruption rather than accretion because of the high 
velocity collision unless their pericenters are aligned (e.g., Safronov 1969; Greenberg et al. 1978; 
Ohtsuki 1993; also see section 2). Then, planetesimal accretion would be forestalled. As shown in 
Eq. (5) in section 2, if pumped-up e and i are larger than 0.01, the velocity dispersion exceeds the 
surface escape velocity of planetesimals in the early stage. Therefore, small orbital modification 
(e, i > 0.01) can give significant influence on planet formation. If the longitudes of pericenter and 
ascending nodes of colliding planetesimals are aligned (" phase alignment"), the colliding velocity 
is significantly reduced from that in the above argument (Marzari and Scholl 2000; also see section 
6). Hence, we also examine the degree of the alignment. 

This work is motivated by Ida, Larwood, and Burkert (2000; hereafter ILBOO). They showed 
that the radial gradient of the pumped-up e and i is rather steep; e and i are highly pumped up in 
the outer planetesimal disk, while the inner disk is almost intact (also see sections 3 and 4). The 
steep gradient leads to a sharp boundary of the disk that divides the strongly perturbed region 
where planet formation is inhibited and the intact region where planet formation keeps going. The 
boundary radius determines radial size of the region of a planetary system where planetary-sized 
bodies exist. Most planetesimals in the outermost region are ejected (and some of them are captured 
by the passing star) , which would truncate the region of the planetary system where solid materials 
exist. Here we are concerned with the former size. 



In the present paper, we investigate the orbital modification due to a passing stellar encounter 
and discuss the effects on planetary formation. We derive analytical formulae of pumped up e 
and i, because, as stated above, the boundary between the planet forming region and the planet- 
formation inhibited region is marked by e, i as small as 0.01 and hence linear perturbation analysis 
is available to determine the boundary. We also did numerical simulations of 10,000 test particles 
in some parameter ranges. Our analytical formulae show excellent agreement with the numerical 
results. 

Encounters between a star and a particle/gas disk have been studied by many authors in 
different contexts. However, most previous studies were concerned with the strongly perturbed 
region in the disk. 

Galactic encounters may form galactic tidal bridges and tails. Many numerical simulations 
(e.g., Toomre and Toomre 1972; Barnes and Hernquist 1992and references therein) have been 
done and analytical calculations with impulse approximation was applied for high speed galactic 
encounters (e.g., Binney and Tremaine 1987). Kalas and Jewitt (1995) suggested that asymmetry in 
the dust disk around /? Pic is caused by a passing stellar encounter. Kalas et al. (2000) and Larwood 
and Kalas (2001) demonstrated through numerical simulation that a passing stellar encounter 
reproduces the observed disk asymmetry and ringlet structure. These structure in galactic and 
dust disks is formed in the outer part of a disk which is strongly perturbed. 

A galaxy may be captured by another galaxy through energy transfer from orbital motions 
to internal motions caused by a close encounter to be a satellite galaxy (Palmer and Papaloizou 
1982; Wahde et al. 1996). Similarly, an encounter of a star with another star with a protoplanetary 
gas disk may lead to capture of the passing star to form binary stars. It may also make a new 
companion from the disk (Boffin et al. 1998), cause disk truncation (Hall 1997), or cause rapid 
disk accretion (Ostriker 1994). The transfer of energy and angular momentum between a passing 
star and the disk is calculated numerically (Clarke and Pringle 1991; Hall et al. 1996; Heller 1993; 
1995) and analytically (Ostriker 1994). Korycansky and Papaloizou (1995) and Larwood (1997) 
did more detailed analytical calculations. These studies were mostly concerned with total energy 
and angular momentum transfer, that is, concerned with strongly perturbed regions in the disk. 

Ostriker (1994) also studied weakly perturbed regime, although the author mainly focused on 
the strongly perturbed region where resonant effects are important. The author only presented 
approximately averaged change in angular momentum in the weakly perturbed region. As we 
discuss later, the author's results would not give enough information to deduce change in e and i, 
which we want to use to discuss the effects on planet accretion. 

Passing stellar encounters are also important for planetary systems or planetesimal disks. The 
planetesimals ejected by jovian planets may become weakly bounded in the Solar system to form 
Oort cloud. Passing stellar encounters would make their binding stronger (Brunini and Fernandez 
1996; Fernandez 1997) or send them back to inner Solar system as long-period comets (Eggers 
and Woolfson 1996; Yabushita et al. 1982). They considered stellar encounters with planetesimals 



ejected from a protoplanetary system. 

At present, mean distance between stars in the solar neighborhood is so large that only Oort 
cloud would be affected by passing stellar encounters. However, since stars are generally born as 
members of an open cluster and stay in the cluster on timescales more than 10^ years (Kroupa 
1995; 1998), passing stellar encounters affect a planetary system or a protoplanetary planetesimal 
system itself in the early stage. Monte-Carlo numerical simulations suggest that the encounters 
can modify nearly circular orbits of giant planets to eccentric orbits, which may correspond to 
observed extrasolar planets in eccentric orbits (de la Fuente Marcos and de la Fuente Marcos 
1997; Laughlin and Adams 1998; 2000; Bonnell et al. 2001). ILBOO showed through numerical 
simulations that a stellar encounter may explain high orbital eccentricities and inclinations of outer 
Kuiper belt objects. They suggested that planetary formation was inhibited by the high orbital 
eccentricities and inclinations in the Kuiper belt while more inner region was almost intact and 
planetary formation proceeded, if the most effective stellar encounter with the proto-Solar system 
had pericenter distance ~ 150-200AU. We perform more extensive calculations by both numerically 
and analytically to discuss the inhibition of planet formation in detail. Heggie and Rasio (1996) 
derived analytical formula for the change in eccentricity of binary stars by an encounter of a single 
passing star, using the Laplace- Runge-Lenz vector. In the limit of null mass of a binary companion, 
their formula is reduced to our analytical formula, although we used different method with the 
Gauss's equations (Brouwer and Clemence 1961). We derive the change in not only eccentricity 
but also inclination and study relative velocity of planetesimals taking into account the alignment 
of longitudes, in order to apply to planetary formation. 

Heppenheimer (1978), Whitmire et al. (1998) and Marzari and Scholl (2000) study planet 
formation in a system with a binary companion orbiting outside a circumstellar disk around the 
primary star. The circumstellar disk is perturbed by the binary companion. They investigated 
relative velocity between planetesimals in the perturbed disk to find the condition of inhibition of 
planet formation, while we consider the effect of one passage of a field star. 

In section 2, we briefly summarize critical velocity dispersions or critical eccentricities of plan- 
etesimals that are important to discuss planet accretion. In sections, we explain calculation models. 
In section 4, we show results of numerical simulations. The analytical formulae are derived in section 
5. e and i pumped up by stellar perturbations in the weakly perturbed inner regime are explicitly 
given as functions of heliocentric radius and parameters of stellar encounters. Using the analytical 
formulae, in section 6, we discuss the size of planet forming region as a function of parameters of 
stellar clusters. We will show the radius of planet forming region is likely to be 40-60AU, in the 
case of a dense cluster, which may be consistent with the Solar system. 



2. PLANETARY ACCRETION 

Planetesimals with 1-10 km sizes are formed from dust grains in a plotoplanetary disk through 
self-gravitational instability (e.g., Safronov 1969; Goldreich and Ward 1973) or sticking induced by 
turbulence of the disk gas (e.g., Weidenschilling and Cuzzi 1993). Accretion of the planetesimals 
make terrestrial planets and cores of jovian planets. The sticking (accretion) of planetesimals 
is caused by self-gravity. The accretion occurs, only when rebounding velocity of planetesimals 
after a collision is smaller than their surface escape velocity (fesc)- Since collision velocity is ~ 
\/Vo+visc (Ohtsuki 1993), where vq is relative velocity of unperturbed crossing orbits neglecting 
mutual gravity, accretion occurs when vq < Vesc; except when too elastic or too inelastic cases 
(Ohtsuki 1993). If stellar encounters pump up vq more than Ueso a collision between planetesimals 
results in disruption rather than accretion (e.g., Safronov 1969; Greenberg et al. 1978). Then 
planet formation is forestalled in this region. If the longitudes of pericenter and ascending node 
are randomly distributed, vq is nearly equal to velocity dispersion given by Wd ~ ve^ + i v^^p 
(e.g. Safronov 1969; Lissauer and Stewart 1993; Ohtsuki et al. 1993). As shown in section 6, the 
longitudes of the perturbed planetesimal orbits are not always random. We thus define "effective" 
eccentricity and inclination, 

Ceflf = bo,a;j/|Mcp, (1) 

^eflf = bo,z|/«kep, (2) 

where uo.xj/ and uo,z are the component on the initial disk and the vertical to the disk, respectively. 
By definition, e^s ^ e and iegf < i. If the longitudes are aligned, icflf = and e^s expresses small 
relative velocity due to shear motion. 

The surface escape velocity of planetesimals with mass m and internal density p is 



1 



(3) 



As shown in section 6, Ccff ^ icflf) so that vq ~ eeflffkep- The condition for a disruptive collision is 
Vq > fcsc = eeff,crit^kep (s-g-. Stern 1995). The condition with effective eccentricity is 



,2^3^\ i 



f 3271171'^ a'^p , ,^, 

Ceflf > eeflf,crit ^ I ^^ I , (4) 

where a and Mi are semimajor axis of a planetesimal and the mass of the host star, respectively. 
Numerically, Ceff^crit is 

,,,„,..„^„l(_^)'(_Z_)^(_^)', ,5) 

where the nominal value of m ~ 10^^ g corresponds to typical mass of Kuiper belt objects. In 
section 6, we show that the phase alignment breaks down for e > 0.01, so that e^s ~ e for e > 0.01 



and the condition (4) is no other than the condition for e. Equation (5) shows that e as smah as 
0.01 can significantly affect planetary accretion. If we consider "protoplanets" with lunar mass, 
Ceff.crit IS ~ 0.1. However, it does not change the discussion on planet accretion in section 7. 



3. CALCULATION MODEL AND BASIC EQUATION 

We model a planetesimal disk as non-self-gravitating, collisionless particles that initially have 
coplanar circular orbits around a primary (host) star, because two-body relaxation time and mean 
collision time between planetesimals are much longer than an effective encounter time scale that 
is comparable to Kepler time scale at pericenter distance (D) of the encounter (for example, it is 
~ 10^ years for D ~ lOOAU). We also neglect hydrodynamical gas drag, because the damping time 
due to the drag [~ W'^ {'m/lO'^^g)^/^{e/0.1)'^^ yr at lOAU (Adachi et al. 1976)] is longer than the 
effective encounter time ~ 10^ years for m > lO^'^g at lOAU; we are considering planetesimals with 
much larger masses. The particulate disk encounters a hypothetical passing star. Note that the 
gravitational relaxation, collision and the drag can be important when we consider planet formation 
on a longer time scale after the stellar passing. The equation of motion of a planetesimal in the 
heliocentric frame (the frame with the primary star at center) is 

fll-_^MLr.+ G^'- (R-r.)-^^R (6) 

dt^ ~ I r, P"^^ + I R - r, |3 ^^ ""^^ I R |3 ^' ^^^ 

where Mi and M2 are masses of the primary and the passing stars, Yj and R are position vectors 
of the planetesimal j and the passing star. The first term in the r. h. s. is force to produce Kepler 
motion around the primary star, and the second and third terms are direct and indirect perturbing 
forces of the passing star. 

We scale length by pericenter distance D of the stellar encounter, mass by the primary star 
mass Ml, and time by ^^cv where ilkcp is Keplerian frequency at a = D given by ^JGMi/D^. 
Equation (6) is then transformed to 



(fvj _ Yj M^{rj-R) M^R 

IF ~ ~\rj |3 " Ir^-R |3 ~ |i?|3' 



(7) 



where M^ = M2/M1, Yj = Yj/D, R = R/-D, and t = ri^cp^- Thus the parameters of encounters are 
inclination (i*) relative to the initial planetesimal disk, eccentricity (e*), and argument of perihelion 
{io^) of orbit of the passing star, and the scaled passing star mass (M*). The encounter geometry 
is illustrated in Fig. 1. We calculate changes in e and i of the planetesimals according to Eq. (6) 
or (7) with various encounter parameters, through orbital integration and analytical estimation. 



4. NUMERICAL SIMULATION 

Regarding the method of numerical integration, we follow ILBOO. We integrated orbits of 
10,000 particles with surface number density Us oc a~^'^. The particles are distributed in the region 
a/D = 0.05-0.8. Since we neglect mutual gravity and collisions of planetesimals, the particular 
choice of a-dependence of Ug and outer and inner edges of the disk does not affect the results. The 
initial e and i of particles are 0. We integrated Eq. (7), using a fourth order predictor-corrector 
scheme. Much more variations of encounter geometry, encounter velocity, and passing star mass 
were examined than ILBOO did. 

Figures 2 show time evolution of e (left panels) and i (middle panels) and corresponding face-on 
snapshots (right panels) in the case with e* = 1 (parabolic orbit), i^ = 30°, w* = 0°, and M^ = 1. 
The (a) top, (b) middle and (c) bottom panels show snapshots at t = —1.33, and 1.33. The 
resultant e and i of planetesimals are mostly acquired when the passing star is near the pericenter. 

The pumped-up e and i are shown in Fig. 3, as a function of the scaled initial semimajor axis 
a/D, in the case with i^ = 5° and uj^ = 90°. (M*, e*) = (1,1), (0.2,1) and (1,5) in Figs. 3a, 3b 
and 3c, respectively. In all cases, we find three characteristic regions of pumped-up e and i. In the 
inner region at a/D < 0.1-0.3, e and i are in proportion to (a/D)^'^ and (a/D)^'"^, respectively. 
In the outer region at a/D > 0.1-0.3, e and i have steeper a-gradient and divergence due to initial 
mean anomaly of planetesimals. In the outermost region, e of many planetesimals, is greater than 
1, that is, the particles are ejected from the system. ILBOO also found these features. We find that 
e oc (a/D)^' ^ and i oc (a/D)^''^ always hold. Analytical estimate derived in section 5 and Appendix 
reproduces e oc (a/D)^''^ and i oc (a/D)'^''^. Note that the dashed lines with triangles are analytical 
estimate derived in section 5 and Appendix. Figures 3 show larger M* and/or smaller e* produce 
larger e and i. The agreement between the numerical and analytical results imply that e and i are 
scaled by M^/ \/M^ + 1. The dependence on e* is a more complicated form (see Appendix, Eqs. 
(46) to (49)). 

Figures 4 show the dependence of i*, in the case with a;* = 0°, e* = 1 and M* = 1. i* are 
5° (Fig. 4a), 30° (Fig. 4b), 45° (Fig. 4c), 85° (Fig. 4d) and 150° (Fig. 4e). An encounter with 
0° < i* < 90° is a prograde encounter relative to rotation of the planetesimal disk and one with 
90° < itf < 180° is retrograde. The numerical results show i is dependent on i^ like oc sin2i^< in 
the inner region. Comparison between Figs. 4b and 4f suggest that e and i are the same between 
itf and 90° — i* in the inner region, although retrograde encounters lead to less steep a/D-gradient 
than that in prograde encounters in the outer region [a/D > 0.1-0.3). This shows that a secular 
effect is at work, so that we take the orbital average of the particle (Larwood 1997) in section 5. In 
Figs. 5 and 6, we show detailed dependence on i^, as well as cj*. In the inner region, e = eo(a/i?)^'^ 
and i = iQ{a/D)^''^. We investigate the dependence of cq and io on uj^, and i^,, in the case with 
e* = 1 and M* = 1. According to the analytical results (Eqs. (15) and (16) in section 5), we scale 
Co and io by (157r/32-v/2)(l/\/2) and (37r/8\/2)(l/\/2), respectively. We denote the scaled eo and 
^0 by eo and io- Figure 5a shows contours of numerically obtained cq as a function of z* and a;*. 



The analytical result (Eq. 15) is plotted in Fig. 5b. We find dependence on uJ.^, is weaker than 
that on i*. Since the analytical result shows that iq is independent of lo^ (Eq. 15), we plot iq 
as a function of i* in Fig. 6. Filled circles in Fig. 6 show numerically obtained iq. We plot the 
numerical results of ten different w* for each i^ in Fig. 6. iq is almost completely independent of 
u;^, which is consistent with the analytical result. Figures 5 and 6 show that the analytical results 
perfectly agree with numerical results. These figures show cq and iq are symmetric with respect to 
i* = 90°. Figure 6 clearly shows i is proportional to sin2i*. 

We next show the range of the power-law inner region is regulated by Jlkcp and fi* (ILBOO 
also discussed this issue), where ri^gp = \JGM\ jo? and fi* = ^G(M\ + M2)(l + e^)/D'^ are the 
Keplerian frequency of a planetesimal with semimajor axis a and the angular velocity of the passing 
star at pericenter. Figures 3 and 4 show that if Jl^cp > ^^* (0/-^^ ~ 0.2 in case with e* = 1 and M* = 
1), e is in proportion to (a/D)'^-^. If Q* < ri^cp ^ 5il* (0.2 < a/D < 0.3 in case with e* = 1 and 
M* = 1), e is pumped up most steeply. In this region, resonant interactions are important (Ostriker 
1994). The condition of n : 1 commensurability is il^/rj^cp = {a/D)^''^y^{l + M*)(l + e*) = 1/n. 
For example, the 5:1, 4:1 and 3:1 resonances are at a/D ~ 0.21, 0.25 and 0.30 in the case with 
e* = 1, M* = 1 and 0° < i* < 90° (the prograde encounters; Figs. 4a to d). The resonances lower 
than 5:1 dominate non-resonant effects in this case. The 1:1, 2:1 and 3:1 are at a/D ~ 0.63, 0.40 
and 0.30 in the retrograde encounter (Fig. 4e). In the retrograde cases, resonances occur with 
particles in the far side of the disk, so that their effects are relatively weak. Numerical simulations 
suggest that the boundary between the inner and outer regions is at the 5:1 commensurability, 
a/D ~ [(1 + M^:){1 + e*)]^^'^(l/5)^'^, in the case of prograde encounters and at the 3:1, a/D ~ 
[(1 + M^,)(l -|- e*)]^^'^(l/3)^'^, in the case of retrograde encounters. Ostriker (1994) analytically 
derived consistent conditions. If f^kcp ^ ^* {o-/D ^ 0.63 in case with e^= = 1 and M^, = 1), there is 
no Lindblad resonance, so that a-gradient is less steep than that in the resonant region. 

In the outer and outermost regions, pumped-up e and i depend not only on initial radial 
position but also on azimuthal position of planetesimals. In the case of a prograde encounter, 
particles in the near side of the disk at pericenter passage are affected by resonances while ones in 
the far side are hardly affected. In these region, r^kcp is not large enough compared with Vt^,, so 
that such asymmetry remains. In the case of a retrograde encounter, only far-side planetesimals 
are affected. Note that particles with very small e and i exist at resonant points depending on the 
azimuthal position, which might be able start runaway accretion (section 6). 

The change in semimajor axis, Aa/a, is much smaller than that in e and i (in radian) in the 
inner region, as shown by comparison of Fig. 3a and Fig. 7. In the inner region, f^kep ^ ^*) 
so that for particles there gravitational potential of the passing star is quasi-stationary. In the 
stationary potential, energy and hence a is conserved. Since the potential is not axisymmetric in 
the heliocentric frame, e is changed. Since it is inclined from the orbital plane of the particles, i 
is also changed. In the outer and outermost regions with il* ~ rikcp, change in Aa/a can not be 
neglected compared with that in e or i 



A stellar encounter can cause disk truncation. Many particle are unbound (e ^ 1) at a/D ^ 0.3 
after the prograde encounter, in the case with M* = 1, e* = 1 (Fig. 2c). Some particles are cap- 
tured by the passing star during a prograde encounter. Their eccentricities are usually close to 
unity in the frame in which the passing star is at center. (If the passing star also has a planetesimal 
disk, some planetesimals are captured by the primary star.) Note that exact amount of the disk 
truncation and the capture depends on the radius of outer edge of the disk relative to D. Such 
strong encounters would have many interesting features as mentioned in Section 1. However, as 
shown in section 2, in order to discuss the effects on planet accretion, the regimes of e, z ~ 0.01 are 
important. The pumped-up e and i in such regimes are near the inner region. Hence, analytical 
linear calculations would well predict the effects on planet accretion. 



5. ANALYTICAL CALCULATION 

We derive analytical formulae of pumped-up e and i in the inner region. Orbital integrations 
show that in that region, e = eo(a/-D)^'^ and i = io(o/-D)^ • The analytical formulae explain 
these dependence on {a/D) as well as dependence of cq and iq on M^,, i*, a;^, and e^,. Heggie and 
Rasio (1996) derived change in eccentricity of binary stars which encounter a passing single star, 
using the Laplace-Runge-Lenz vector. Using the Gauss's equations, we derive change in inclination 
as well as eccentricity. In the limit of null mass of a binary companion, the formula by Heggie 
and Rasio (1996) is equivalent to out formula for eccentricity. The analytical results show the 
phase alignment, so that e and i do not necessarily determine relative velocity between perturbed 
planetesimals. We numerically study the degree of the phase alignment in section 6. 

We are mostly concerned with parabolic encounters {e^, = 1), since for encounters we are 
interested in, e^, is not far from 1, as follows. (We also calculate hyperbolic encounters (e^, > l)with 
the same analytic method in Appendix.) The specific angular momentum (/*) and specific energy 
{E^) of the passing star orbit relative to the primary star are ~ f*D ~ A/^*d + [^kcp(-D)]^ ^ 
d/ 



D and f*j/2, respectively, where v*^ is velocity dispersion of stars in a cluster. Thereby, e^, 
Vl + 2/2E,/G2(Mi + M2)2 ~ v^l + (t;*dMep)2[l + (t;*dMep)]2. Since v*^ is ~ 1 km/s (Binney 
and Tremaine 1987) and fkep(-C^) is larger than 1 km/s for D < lOOOAU, which we are interested 
in, e* is ~ 1. 

We adopt the following approximations to derive pumped-up e and i in the inner region. 

(i) Aa/a is neglected, since Aa/a <^ e, i, 

(ii) orbital averaging is applied for planetesimal orbits, since ilkcp S> ^2* (Larwood 1997), 
(iii) e, i, a/D <^ 1. 

In the equation of motion of a planetesimal given by Eq. (6), the second and third terms in 



10 



the r. h. s. stand for the perturbation forces of the passing star. We define 

GM2 ,^ GMa^ 

Fperturb - 1 p^ _ ^ 13 (R " r) - Tp-pl^- (8) 

We divide the perturbation force into r, 9, and z components, 

Vpcrturh = R^r + feQ + Ne^, (9) 

where e^ (= r/ | r |), e^ and e^ are unit vectors in the radial, tangential, and normal components 
in initial orbital plane of the planetesimal disk, respectively. We define h = esin(ti; + ri),A; = 
ecos(tiJ + Vt),p = sinisinfi and q = sin i cos Q, where uj and Vl are the argument of pericenter and 
longitude of ascending node of a planetesimal. The equations of motion with these orbital elements 
in h,k,p,q <^ 1 {e,i <^ 1: assumption (iii)), which are called Gauss's equations, are (Brouwer and 
Clemence 1961) 



dh I a (_jjcos0 + 2rsin^), (10) 



dt V GMi 

dk 



dt V GMi 



"' {Rsme + 2Tcos9), (11) 



dp _, / « -iVsin^, (12) 



dt V GMi 



'^'^ - ^^~-Ncose, (13) 



dt V GMi 

where 9 = f + uj + 0, (/is true anomaly) and we retained the lowest order terms of e and i in the 
right hand side. 

We expand R, T, N in terms of a/D up to the terms of {a/D)"^ (see Appendix, Eqs.(37), (38), 
and (39)). Corresponding to assumption (ii), we take orbital averaging, e.g., 

•dh\_j^''{dh/dt)d9_ ^^^^ 



dt / 27r 

The averaged dh/dt, dk/dt, dp/dt and dq/dt are shown in Eqs. (41), (42), (43) and (44) in 
Appendix. Finally, we integrate them with the stellar passage, assuming constant a of planetesimals 
(assumption (i)) to obtain changes in h, k, p, and q caused by the stellar passage in parabolic orbit 
(Appendix, Eqs. (46) to (49)). 

Since we start with h, k,p,q = 0, their changes are equal to final h, k, p, and q of planetesimals. 
Since e = yh?-\-Iv and i = \Jp^ + (7^, we obtain final e and i as 

ISvr M. 



32^/2 ^/TTM* ^D' 

cos^ cj* I 1 - - sin^ i* I + sin^ a;* cos^ i* f 1 — sin^ i* ) , (15) 
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37r 



M, 



sin2i*|, 



8^2 vttm; VL>. 

cosaj*(5cos^ i* — 1) 



arctan 



n 



sinw^, cosi^,(15cos2 i^ — 11) 

7r/2 for 0° Si*S 90°, 

-7r/2 for 



90° <uS 180°, 



(16) 

(17) 

(18) 



where oj is the longitude of pericenter and = lo + 0,. The detailed derivation is shown in Appendix. 
We have analytically explained the power-law dependence of pumped-up e and i on (a/D) in inner 
region. Since the orbital averaging is valid, Co and O of the perturbed orbit do not depend on 
azimuthal position of a planetesimal as Eqs. (17) and (18). Such phase alignment of planetesimals' 
orbits, leads to reduction of the colliding velocity between planetesimals (section 6). 

We plot the analytical formulae (15) and (16) with dashed lines with triangles in Figs. 3 and 
4. eo (= e/{a/ D)^'"^) and i^ (= i/{a/D)^''^) of the analytical formulae are showed in Figs. 5b 
and 6. The analytical formulae show excellent agreement with the numerical results in the inner 
power-law inner region. We examined the agreement in the cases with other parameters and found 
that the difference is always less than a factor 2. 

The pumped-up e and i depend on phase angles i^, and oj* of a stellar encounter. The average 
of e and i with cj* and i* is 



1 

47r 

" 1 

47r 



2tt 



dujt 



di^(s'mi^)e 



2tt 



cLl)^ I (ii*(sini*)i^ 



ISvr 



M* 



32^7 yiTM* 
\/l57r M^ 



20 VM, + 1 
We use Cave and iave to estimate the boundary of planet accretion in section 7. 



(19) 
(20) 



Ostriker (1994) analytically investigated the change in energy (AE'disk) and perpendicular an- 
gular momentum (ALjisk^^) of an entire circumstellar disk during a parabolic stellar encounter. 
Although the author was concerned mainly with resonant region for rdiskmax ^ 0.2Z), the au- 
thor presented phase-averaged approximate ALdisk,^ (but not AE'jisk) for rdisk.max ^ 0.2D (non- 
resonant region). Since specific perpendicular angular momentum AL^ are written by Lo(Aa/a 
— Ae^ — Ai^)/2 where Lq is the initial angular momentum, ALdisk,^ is obtained by integrating AL^ 
over radius from rdisk.min to rdisk,max, using our results in the power-law inner region. Our ALdisk,2 
is larger by a factor ~ 0(10) than the author's with less steep a-gradient. 

Heggie and Rasio (1996) derived formula of change in eccentricity of binary stars after a single 
passing encounter. To derive this formula, they use the Laplace-Runge-Lenz vector. Their formula 
of change in eccentricity of binary stars with circular orbit is perfectly the same as our formula 
derived by a different method, in limit of null mass of a binary companion. 
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6. RELATIVE VELOCITY 

As shown in the last section, difference of uj or 0, between planetesimals are almost as long 
as the linear analysis is valid. We investigate the relative velocity between planetesimals after a 
stellar passing, using the numerical data. Figure 8 show that a) and 0, after the stellar encounter 
with (i*, uj^:, M*, e*) = (5°, 90°, 1, 1). lu and fi are aligned in the region where e and i have 
the power law radial dependence. In the outer region where the resonances are important, u) is 
randomly distributed (Retrograde encounters show similar results). 

Whitmire et al. (1998) evaluated relative velocities between planetesimals whose orbits are 
crossing assuming coplanar orbits. Extending their method to three-dimension with assumption 
i < l(radian), we evaluated the relative velocity between planetesimals 1 and 2. For « ^ 1, the 
projective orbit of a planetesimal on the initial disk is almost the same as the case with i = 0. We 
calculate the relative velocity, if \zi — Z2I is smaller than the sum of planetesimal radii at the point 
where the projective orbits are crossing, where Zi(2) is position vertical to the initial disk. The 
component on the initial disk VQ^xy and that vertical to the disk vq^z at a crossing point are 



'^0,xy 



VO,; 



GM^ 



ei sin(/2 - AlD) 62 sin /2 



VGMi 



Pi 
h[cos{f2 



+ 



Aoj + uJi) + eicoswi] i2[cos(/2 + UJ2) + 62 cos ^2] 



(21) 
(22) 



where pi(2) = Oi(2)(l~ei(2))) AcD and AQ are u)2—^i and ^22 — ^^1. /2 is true anomaly of planetesimal 
2 at the crossing point, satisfying 



cos/2 



-AB ± CVC^ + B^-A^ 



(23) 



where 



A = P2-P1, 

B = eip2 cos Auj — e2Pi, 

C = eip2sinAa). 



(24) 
(25) 
(26) 



r is radial position of planetesimals at the crossing point, given by 

P2 

r = . 

1 + 62 cos /2 



(27) 



In Fig. 9, e and i are compared with Ceff = I'o.xy/'Wkep and igflf = ^o,^/^kcp) which are calculated 
through Eqs. (21) and (22) with the data in Fig. 8. We also plot 



Brci = |e2 — ei| = e^ + 62 — 2eie2 cos Aw, 
^rci = lb — ill = ^1 + ^2 ~ 2iii2 cos A$7, 



(28) 
(29) 
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where e^ ij (j = 1,2) are (cj cos ujj , Cj sin Coj ) and {ij cos Vtj , ij sin 0^ ) . In Fig. 9, we found ej-ci — eeS 
and ij-ci — ^cflf- The effect of the phase ahgnment is more clear with Crei and irei rather than with 
Ceff and Zeff, because e^eX and irei are directly related with Aci) and Afi. A stellar encounter results 
in alignment of Cj and Jl (Aa) ~ and AO ~ 0) in inner region, e^s (crei) "^ e and iefr (irei) "^ i- In 
outer region, the linear approximation breaks down and Au; and AJ7 are randomly distributed, so 
that Crci — e and irci — i- Figure 8 shows that the alignment breaks down at a/D ~ 0.2 for uj, and 
~ 0.3 for $7, resulting in more rapid increase in Crci with a than ij-ci- As a result, at Ccs ~ ecff,crit ~ 
0.01, Bcff ~ e > ics- Also for other encounter parameters, we find a similar trend. We thus estimate 
planet formation region, using Eqs. (5) and (15) with e^s ~ e. 

In Fig. 9, we also find that at a resonant point {a/D ~ 0.2), all of e and e^s have wide variety of 
values caused by the resonant feature. Although many planetesimals have Ces well beyond ecff,crit, 
but some fraction has Cefr <C eeflf,crit should be very fast. If it can precede disruption due to collisions 
with high Ccff planetesimals to form large bodies with high eefr,crit quickly, planet can be formed at 
this particular resonant point. 

We also study the cases with non-zero initial e and i because self-gravity of planetesimals or 
more distant stellar encounters would occur before the strongest encounter. Figure 10 shows the 
e, i, oj and il of planetesimals after the encounter with initial e and i = 1 x 10~^. lD and Q are 
initially at random. The encounter parameters are (i^,, a;*, M*, e*) = (5°, 90°, 1,1). If change of e 
and i by the encounter are much larger than the initial e and i, the features of e and i are almost 
the same as those in the case of initial e and i = 0; oj and fi are aligned If the pumped-up e and i 
are smaller than the initial e and i, uj and fi remain at random. 

After the stellar encounter, self-gravity of planetesimals, collision between them and gas drag 
would be important for accretion of planetesimals on a longer time scale. They may keep the phase 
alignment rather than destroy it (e.g., Goldreich and Tremaine 1982; Marzari and Scholl 2000; Ito 
and Tanikawa 2001), until a strong perturber such as a giant planet is formed. We need more study 
about these effects on planetesimal accretion. 



7. SIZE OF A PLANET FORMING REGION 

In the last section, we have derived analytical expressions for e and i pumped up by a passing 
stellar encounter. The expressions perfectly agree with numerical calculations in the region e < 0.01. 
In this section, we estimate the size of planet forming region in a protoplanetary disk, using these 
analytical expressions. 

According to the argument in the last section, the collision between planetesimals results in 
disruption rather than accretion (Eq. (5)), 



«>^~°-"'llS^J imuJ • <^°> 
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where m is mass of a planetesimal. For e > 0.01, oj and fi are not aligned and e ~ e^s- 

As shown in sections 3 and 4, the radial gradient of e is so steep that there is a sharp boundary 
of the disk that divides the planet forming region and the disruptive region where planet formation 
is inhibited. Since Ces, which includes the effect of the phase alignment, has steeper gradient than 
e, the actual boundary should be sharper. Here we derive the radius of the boundary as a function 
of physical parameters of the stellar cluster which the host star of the planetary system belonged 
to. 

In the case with e > 0.01, the analytical expressions slightly underestimate the pumped-up 
e, since resonant effects are also important in this case. However, we can use them to estimate 
the boundary radius, because the radial gradient of e is so steep that the underestimation hardly 
changes the boundary radius. 

In a stellar cluster, e^ ~ 1, as mentioned in Section 5. Substituting Eq. (19) into Eq. (5), we 
obtain the boundary radius Cpianet as 



m 
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6 / F\ 4 / D \ 4 



( 



<'*»~''%I0^J [jj 1,I50AUJ ^"- <=*'> 

where F = M* + 1/M^. If M* is 1, F is 2. The factor (F/2)^/'^ can not significantly deviate from 
1. The dependence of m on apianct is also very weak. Therefore, apianet depends almost only on D. 
A stellar encounter with D ~ 150-200AU restricts the disk radius of a planetary system (the disk 
radius of planet forming region) to 40-60AU. 

Considering evaporation process of a stellar cluster, Adams and Laughlin (2001) estimated 
effective D before the evaporation as 

D ~ 200 f :5£!H!!££^ ( JL] aU, (32) 

V 2pc J V2oooy ^ ^ 

where N is number of stars in a stellar cluster and -Rciustcr is size of the cluster. For the Trapezium 
cluster in Orion, N ~ 2300 and -Rciustcr ~ 2 pc. In a dense cluster like Orion Trapezium, D is 
as small as 200 AU, so that Opianet ~ 40-60 AU. ILBOO demonstrated that the high eccentricity and 
inclination of objects in the outer Kuiper Belt may be explained by the stellar encounter with 
D ~ 150-200AU, which may suggest the Sun was born in a dense stellar cluster. 

So far, we have only considered passing of a single star, however, passing of binary stars would 
also be important. Laughlin and Adams (1998) and Adams and Laughlin (2001) (also see the 
next section) suggested that passing binary encounters are more disruptive than passing single-star 
encounters. If we take into account the effects of passing binary encounters, apianct rnay be smaller 
than Eq. (31). 

We should investigate distribution of D (not only an effective value) as well as the effects of 
passing binary encounters, to discuss diversity of sizes of planetary systems in more detail. 
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Note that as discussed in the last section, less than planetesinials might coalesce with each 
other, resulting in a kind of runaway growth; planets might be formed at particular resonant 
location beyond apianet- 



8. CONCLUSION & DISCUSSION 

We have investigated the effects of a passing stellar encounter on a planetesimal disk through 
orbital integration and analytical calculations. Since stars are generally born as members of a 
cluster and would stay in the cluster on timescales more than 10^ years (Kroupa 1995; 1998), a 
relatively close encounter, e.g., one with distance ~ 200AU is likely to occur during formation age 
of a planetary system. 

We considered that a disk of massless particles (planetesimals) orbiting a primary star en- 
counters a passing single star. Encounter parameters are pericenter distance of the encounter {D), 
the argument of perihelion (lu^), eccentricity (e^,) and inclination (i*) of the orbit of the passing 
star, and the mass ratio (M*) of the passing star's mass to the primary one. We showed that the 
pumped- up orbital eccentricities e and inclinations i of planetesimals have steep radial gradient. 
In the inner region at semimajor axis a ^ a [(1 + M^,)/2]^'^ [(1 -|- e^)/2]^'^ D, (a ~ 0.2 for a 
prograde encounter and a ~ 0.3 for a retrograde encounter) e oc {a/ D)^''^ and i oc (a/D)^'^, inde- 
pendent of the encounter parameters. The result is also independent of initial azimuthal position 
of planetesimals, because orbital period of planetesimals is much shorter than passing time scale. 

In the outer region a > a [(l-|-M*)/2]^/^ [(l-|-e*)/2]-'^'^ D, the radial gradient is steeper, but is 
not expressed by a single power-law. In this region, resonant effects are more important (Ostriker 
1994). 

The longitude of pericenter Co and ascending node Vt are aligned between planetesimals, in 
power-law region of e and i (Fig. 8). We investigated relative velocity between planetesimals whose 
orbits cross. We define CcS = 'Wo,a;j//^kcp and i^s = vo^z/vkep where vo^xy and vq^z are planar and 
vertical components of the relative velocity. We found that Ces and ies are smaller than e and i by 
many orders of magnitude. In outer region, the phase alignment brake down, so that e ~ Ces and 

i ~ ieS- 

The stellar perturbations significantly affect the outer part of a planetesimal disk. Even in 
the inner disk, Ccs and ics that are pumped up to > 0.01 inhibit further planetesimal accretion, 
because the relative velocity corresponding to Ccg, icS ^ 0.01 exceeds the surface escape velocity 
of planetesimals and collisions between planetesimals should be disruptive. We investigate the 
relatively weakly perturbed region with e, z ~ 0.01 to study the effects on planetesimal accretion. 
We have derived analytical expressions of e and i in the power-law inner region. 

With the approximations of constant the semimajor axis a and orbit averaging of planetesimals' 
motion, we derived analytical expressions of the pumped-up e and i (Eqs. (15) and (16)) in the 
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case with a parabolic encounter (e* = 1) as 



2 /i 5-2-\,-2 2 ■ ( -, 15. 2- 

COS o;* I 1 — — sm «* ) + sm lj^< cos i* | 1 ^ sm i.. 



1/2 

(33) 



? — — 1= , =\-r] sin2z* . (34) 

Note that the dependence on lo^ is weak in both expressions and that on i^ is also weak in the 
former expression (Figs. 5 and 6). The analytical formulae agree with the numerical results within 
a factor 2. We also calculated e and i in the case with a hyperbolic encounter (e^, > 1) in Appendix. 

At Ccff ~ 0.01, e ~ Ceff and Cefr is larger than iefj, because the phase alignment for Co breaks 
down at smaller a than that for i7. Hence the relative velocity is regulated by e at e^s ~ 0.01 and 
disruptive collisions e > eeflf,crit where eefr,crit is given as Eq. (5). Since the radial gradient of e and 
Ceff are so steep that there is a sharp boundary in the disk that divides the planet forming region 
and the disruptive region where planet formation is inhibited. Planetesimal orbits are significantly 
modified beyond the boundary, while they are almost intact inside the boundary. We derived the 
boundary radius (apianct) by the condition that the pumped- up velocity dispersion of planetesimals 
is equal to their surface escape velocity (Eq.(31)): 

„„,„-.40(^)"(f)'(^)'AU, (35) 

where F ~ M* + 1/M^ and m is planetesimal mass. The variation of {F/2Y''^ is small. The 
dependence of apianct on ra is very weak. Thereby flpianct depends almost only on D. In a dense 
cluster like Orion Trapezium, D is as small as 200AU (Adams and Laughlin 2001). Such a stellar 
encounter with 150-200AU restricts the disk radius of a planetary system (the disk radius of planet 
forming region) to 40-60AU. 

In the Solar system, there would be no planetary-sized object beyond Neptune at 30AU. 
Kuiper-belt objects beyond Neptune have velocity dispersion considerably larger than their surface 
escape velocity. These might be accounted for by a stellar encounter in a dense stellar (ILBOO; 
Adams & Laughlin 2001). Our Solar system may have belonged to a dense stellar in the formation 
age. 

In a dense stellar cluster, planetary systems cannot be significantly larger than the size of the 
planetary region of our Solar system (~ 30-40AU) . To discuss diversity of sizes of planetary systems 
in detail, we will need to investigate distribution of D (not only an effective value). In the present 
paper, we have only considered passing of a single star, however, passing of binary stars would 
also be important. Laughlin and Adams (1998; 2000) suggested that passing binary encounters are 
more disruptive than passing single-star encounters. We performed several runs of passing close 
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binary stars. If the total binary mass is equal to a single passing star mass, the results show that 
the encounter of close binary stars is similar to the single stellar encounter, except that e is pumped 
up more highly at some resonant positions corresponding to the binary frequency even in the inner 
region. As long as close binary cases are considered, apianct may be estimated with the single star's 
equation using binary total mass. 

The effects of cumulative distant (large D) encounters may be neglected compared with a 
few closest encounters. The successive distant encounters change h,k,p,q defined in section 5 
constructively or destructively, so that their averages would be zero. Their dispersions depend on 
D with large negative power-indexes (< —3), so that the cumulative effects integrated by 2iTDdD 
would quickly vanish with D. 

Circumstellar dust disks around stars are observed. This dust may be secondary, because 
dust would be removed before stellar age, by radiative pressure and pointing Robertson effect 
(e.g., Artymowicz 1997; Backman and Paresce 1993). Disruptive collisions at a > apianot would 
continuously produce dust materials, which might form dust-debris disks around Vega-like stars. 
On the other hand, planetesimals grow to be large bodies at a < Opianet; not producing dust 
materials. Although radiative pressure and pointing Robertson effect would cause migration of the 
produced dust materials (e.g., Artymowicz 1997; Backman and Paresce 1993), Opianct correspond 
to the radius of inner holes of dust-debris disks. The radii of the observed inner holes around e 
Eridani (Greaves et al. 1998), HD 14156 (Andrillat et al. 1990; Sahu et al. 1998; Augereau et al. 
1999; Weinberger et al. 1999; 2000), HD 207129 (Jourdain de Muizon et al. 1999) and HR 4796A 
(Jura et al. 1995; Fajardo-Acosta et al. 1998; Jayawardhana et al. 1998) are 30-lOOAU, which are 
comparable to flpianct m the case of relatively dense clusters. We need more detailed analysis of 
distribution of Cpianet as a function of environment parameters of stellar clusters, combined with 
discussions of collision outcomes and migration processes of dust materials, to discuss the diversity 
of inner holes of dust-debris disks. 



Appendix 

We integrate Eqs. (10), (11), (12) and (13) and derive expressions of e and i pumped up by a 
passing star as functions of a/D, M*, lo^, e* and i*. 

Fperturb defined by Eq. (8) is expanded as a series of a/D up to the order (a/D)"^: 



' perturb — pj 



R-e,, \/a\ r/15(R-e^)2 3\ , , ]/a\' 



(36) 

where -R =| R | and e^ = r/a. We denote Fperturb = R^r +Teg +Nez, where eg and e^ are 
unit vectors in the tangential and vertical directions. We denote that e^ = (cos 0, sin ^,0) and 
R = {Rx,Ry,Rz) in the Cartesian coordinates where the initial planetesimal disk is on the x-y 
plane and the x-axis is directed to ascending node of the passing star's orbit. The components. 
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(39) 



where (p \s z component of e^ x R. 

Substituting Eqs. (37), (39), and (38) into Eqs. (10), (11), (12) and (13) and taking the orbit 
averaging, e.g.. 



dh\ _ fQ{dh/dt)de 
It/ ~ 2^r ' 



where "< >" means the orbit averaging, we obtain 
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(40) 

(41) 

(42) 

(43) 
(44) 



Finally, we integrate them along the trajectory of the passing star: 



Rx 

Ry 

R. 



i?cos(/* + w*) 
i?sin(/* +a;^,)cosi^, 
i?sin(/* + LO^) sini* 



(45) 



where f^, is the true anomaly of the passing star. We integrate Eqs. (10) to (13) over time t from 

—00 to 00. 

We integrate (dh/dt), (dk/dt), {dp/dt) and (dq/dt) by f^, instead of by t, using df^f/dt = 

--df* 
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19 



Ap 



M* 



8(e, + 1)3/ VM, + 1 \DJ 



a\3/2 

— (27 + /i) sin 2i, 



e* /l-e-2\^/2 M, /a\3/2 . , . ^ 
Aq ~ — - -— sm i* sm ioj* 

where 7 is cos^^(— 1/e^,) and 

^2 



/ = V^ 



1 



9 
h 



VW^ 



^(1 + 2e2) - (2 + el + Uet) sin^ i, + -(e^ - 1)^ cos 2tu, sin^ i, 

D 24 D 

^(1 + 2e2) + ;^(2 - I9el - 28e^) sin^ i, + hel - if cos 2u;, sin^ i, 
D 24 



-e*Yl-e* {(e^ - 1) cos2l<j* + 3}. 



(48) 
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(50) 

, (51) 
(52) 



Since we start with h, k,p,q = 0, the changes are equal to final h, k, p and q of the planetesimals. 
We obtain final e and i with e = \/h? + fc^ and i = \Jjp' + cp'. For e^, = l, f = g = h = Aq = 0, e 
and i are 

157r M* / a \ 2 

e ~ 



(53) 
(54) 



128^/2 v/TTM* VD 



•*{5cos i* — 1} + sin cj* cos i*{l5cos i* — ll}" 



cos CO: 



3a M* /a\l . ^. 
— sm 2z=. 



8v/2 VI + M* ^^ 



ACKNOWLEDGMENTS 

We thank Kiyoshi Nakazawa and Hiroyuki Emori for useful scientific discussions. We also 
thank John D. Larwood and Giovanni B. Valsecchi for critical comments, which helped us improve 
the original version of the manuscript, and Andeas Burkert for motivating us to study this problem. 
This work is supported in part by Grant-in- Aid for Specially Promoted Research of the Japanese 
Ministry of Education, Science, and Culture (11101002). 



REFERENCES 

Adachi, I., Hayashi, C., and Nakazawa, K. 1976. The gas drag effect on the elliptical motion of a 
solid body in the primordial solar nebula. Prog. Theor. Phys. 56 1756-1771. 

Adams, F. C., and G. Laughlin 2001. Constraints on the Birth Aggregate of the Solar System. 
Icarus 150 151-162 

Andrillat, Y., M. Jaschek, and C. Jaschek 1990. The infrared spectrum of HD 141569. Astron. 
Astrophys. 233, 474-476. 



20 



Artymowicz, P. 1997. Beta Pictoris: an Early Solar System? Annu. Rev. Earth Planet. Sci. 25, 
175-219. 

Augereau, J. C, A. M. Lagrange, D. Mouillet, and F. Mnard 1999. ST/NICM0S2 observations of 
the HD 141569 A circumstellar disk. Astron. Astrophys. 350, L51-L54. 

Backman, D. E., and F. Paresce 1993. Main-sequence stars with circumstellar solid material: the 
Vega phenomenon. In Protostars and planets III (E. H. Levy and J. Lunine, Eds.), pp. 
1253-1304. Univ. of Arizona Press, Tucson. Astrophys. J. 450, 1253-1304. 

Barnes, J. E., and L. Hernquist 1992. Dynamics of interacting galaxies. Annu. Rev. Astron. Astro- 
phys. 30, 705-742. 

Beckwith,S. V. W., and A. I. Sargent 1996. Circumstellar disks and the search for neighbouring 
planetary systems. Nature 383, 139-144. 

Binney, J., and S. Tremaine 1987. Galactic Dynamics. Princeton Univ. Press, Princeton. 

Boffin, H. M. J., S. J. Watkins, A. S. Bhattal, N. Francis, and A. P. Whitworth, 1998. Numerical 
simulations of protostellar encounters - I. Star-disc encounters. Mon. Not. R. Astron. Soc. 
300, 1189-1204. 

Bodenheimer, P., and J. B. Pollack 1986. Calculations of the accretion and evolution of giant planets 
The effects of solid cores. Icarus, 67, 391-408. 

Planetary dynamics in stellar clusters. Mon. Not. R. Astron. Soc. 322, 859-865. 

Brouwer, D., and G. M. Clemence, 1961. Methods of Celestial Mechanics. Academic Press, New 
York. 

Brunini, A., and A. Fernandez, 1996. Perturbations on an extended Kuiper disk caused by passing 
stars and giant molecular clouds. Astron. Astrophys. 308, 988-994. 

Clarke, C. J., and J. E. Pringle 1991. Star-disc interactions and binary star formation. Mon. Not. 
R. Astron. Soc. 249, 584-587. 

Clarke, C. J., and J. E. Pringle, 1993. Accretion disc response to a stellar fly-by. Mon. Not. R. 
Astron. Soc. 261, 190-202. 

de la Fuente Marcos, C. and R. de la Fuente Marcos, 1997. Eccentric giant planets in open star 
clusters. Astron. Astrophys. 326, L21-L24. 

Eggers, S., & M. M. Woolfson 1996. Stellar perturbations of inner core comets and the impulse 
approximation. Mon. Not. R. Astron. Soc. 282, 13-18. 

Fajardo-Acosta, S. B., C. M. Telesco, C. M., and R. F. Knacke 1998. Infrared Photometry of beta 
Pictoris Type Systems. Astron. J. 115, 2101-2121. 



21 



Fernandez, A. J. 1997, The Formation of the Oort Cloud and the Primitive Galactic Environment. 
Icarus 129, 106-119. 

Goldreich, P., and S. Tremaine 1982. The dynamics of planetary rings. Annu. Rev. Astron. Astro- 
phys. 20 249-283. 

Goldreich, P., and W. R. Ward 1973. The Formation of Planetesimals. Astrophys. J. 183, 1051-1062 

Greaves, J. S., and 10 colleagues, 1998. A Dust Ring around epsilon Eridani: Analog to the Young 
Solar System. Astrophys. J. 506, L133-L137. 

Greenberg, R., W. K. Hartmann, C. R. Chapman, and J. F. Wacker, 1978. Planetesimals to planets 
- Numerical simulation of collisional evolution. Icarus 35, 1-26. 

Hall, S. M. 1997. Circumstellar disc density profiles: a dynamic approach. Mon. Not. R. Astron. 
Soc. 287, 148-154. 

Hall, S. M., C. J. Clarke, and J. E. Pringle 1996. Energetics of star-disc encounters in the non-linear 
regime. Mon. Not. R. Astron. Soc. 278, 303-320. 

Hayashi, C. 1981. Structure of the solar nebula, growth and decay of magnetic fields and effects of 
magnetic and turbulent viscosities on the nebula. Prog. Theor. Phys. Suppl. 70, 35-53. 

Hayashi, C, K. Nakazawa, and Y. Nakagawa 1985. Formation of the solar system. In Protostars 
and Planets II (D. C. Black and M. S. Matthews, Eds.), pp. 1100-1153. Univ. of Arizona 
Press, Tucson. 

Heggie, D. C, and F. A. Rasio 1996. The Effect of Encounters on the Eccentricity of Binaries in 
Clusters. Mon. Not. R. Astron. Soc, 282, 1064-1084. 

Heller, C. H. 1993. Encounters with protostellar disks. I - Disk tilt and the nonzero solar obliquity. 
Astrophys. J. 408, 337-346. 

Heller, C. H. 1995. Encounters with Protostellar Disks. II. Disruption and Binary Formation. As- 
trophys. J. 455, 252-259 

Heppenheimer, T. A. 1978. On the formation of planets in binary star systems. Astron. Astrophys. 
65, 421-426. 

Ida, S., J. Larwood, and A. Burkert 2000. Evidence for Early Stellar Encounters in the Orbital 
Distribution of Edgeworth-Kuiper Belt Objects. Astrophys. J. 528, 351-356. 

Ikoma, M., K. Nakazawa, and H. Emori 2000. Formation of Giant Planets: Dependences on Core 
Accretion Rate and Grain Opacity. Astrophys. J. 537, 1013-1025. 

Ito, T., and Tanikawa, K. 2001. Stability of Terrestrial Protoplanet Systems and Alignment of 
Orbital Elements. P^^J 53, 143-151. 



22 



Jayawardhana, R., S. Fisher, L. Hartmann, C. Telesco, R. Pina, and G. Fazio 1998. A Dust Disk 
Surrounding the Young A Star HR 4796A. Astrophys. J. 503, L79-L82. 

Jura, M., A. M. Ghez,R. J. White, D. W. McCarthy, R. C. Smith, and R G. Martin 1995. The fate 
of the sohd matter orbiting HR 4796A. Astrophys. J. 455, 451-456. 

Jourdain de Muizon, M., and 10 colleagues, 1999. A very cold disc of dust around the GOV star 
HD 207129. Astron. Astrophys. 350, 875-882. 

Kalas, P., J. Larwood, B. A. Smith, and A. Schultz 2000. Rings in the Planetesimal Disk of /? 
Pictoris. Astrophys. J. 530, 133-137. 

Kalas, P., and D. Jewitt 1995. Asymmetries in the Beta Pictoris Dust Disk Astron. J. 110, 794-804. 

Korycansky, D. G., and J. C. B. Papaloizou 1995. The response of a gaseous disc to a binary 
encounter Mon. Not. R. Astron. Soc. 274, 85-98. 

Kroupa, P. 1995. The dynamical properties of stellar systems in the Galactic disc. Mon. Not. R. 
Astron. Soc. 277 1507-1521. 

Kroupa, P. 1998. On the binary properties and the spatial and kinematical distribution of young 
stars. Mon. Not. R. Astron. Soc. 298, 231-242. 

Larwood, J. D., and P. G. Kalas 2001. Close stellar encounters with planetesimal disks: The 
dynamics of asymmetry in the f3 Pictoris system. Mon. Not. R. Astron. Soc. 323, 402-416. 

Larwood, J. D. 1997. The tidal disruption of protoplanetary accretion discs. Mon. Not. R. Astron. 
Soc. 290, 490-504. 

Laughlin G., and F. C. Adams 1998. The Modification of Planetary Orbits in Dense Open Clusters. 
Astrophys. J. 508, L171-L174. 

Laughlin G., and F. C. Adams 2000. The frozen Earth: Binary scattering events and the fate of 
the Solar. Icarus, 145, 614-627. 

Lissauer, J. J., and G. R. Stewart, 1993. Growth of planets from planetesimals. In Protostars and 
planets III (E. H. Levy and J. Lunine, Eds.), pp. 1061-1088. Univ. of Arizona Press, Tucson. 

Marzari, F., and H. Scholl 2000. Planetesimal Accretion in Binary Star Systems. Astrophys. J., 
543, 328-339. 

Mizuno, H. 1980. Formation of the Giant Planets. Prog. Theor. Phys., 64, 544-557. 

Ohtsuki, K., S. Ida, Y. Nakagawa, and K. Nakazawa 1993. Planetary accretion in the solar gravi- 
tational field. In Protostars and Planets III (E. H. Levy & J. Lunine Eds.), pp. 1089-1107. 
Univ. Arizona Press, Tucson. 



23 



Ohtsuki, K. 1993. Capture probability of colliding planetesimals - Dynamical constraints on accre- 
tion of planets, satellites, and ring particles. Icarus 106, 228-246. 

Ostriker, E. C. 1994. Capture and induced disk accretion in young star encounters. Astrophys. J. 
424, 292-318. 

Palmer, P. L. & J. Papaloizou, 1982. Tidal interactions of disc galaxies. Mon. Not. R. Astron. Soc. 
199, 869-882. 

Safronov, V. S. 1969. Evolution of the Protoplanetary Cloud and Formation of the Earth and 
Planets. Nauka Press, Moscow. 

Sahu, M. S., J. C. Blades, L. He, D. Hartmann, M. J. Barlow, and I. A. Crawford 1998. Atomic 
and Molecular Interstellar Absorption Lines toward the High Galactic Latitude Stars HD 
141569 and HD 157841 at Ultra-High Resolution. Astrophys. J. 504, 522-532. 

Stern, S. A. 1995. Collisional time scales in the Kuiper disk and their implications. Astrophys. J. 
110, 856-868. 

Toomre, A., and J. Toomre 1972. Galactic Bridges and Tails. Astrophys. J. 178, 623-666. 

Wahde, M., K. J. Donner, and B. Sundelius 1996. Dynamical friction in disc galaxies with non-zero 
velocity dispersion. Mon. Not. R. Astron. Soc. 281, 1165-1182. 

Weidenschilling, S. J., and J. N. Cuzzi 1993. Formation of planetesimals in the solar nebula. In 
Protostars and Planets III (E. H. Levy and J. I. Lunine Eds.), pp. 1031-1060. Univ. of 
Arizona Press, Tucson. 

Weinberger, A. J., E. E. Becklin, G. Schneider, B. A. Smith, P. J. Lowrance, M. D. Silverstone, 
B. Zuckerman, and R. J. Terrile 1999. The Circumstellar Disk of HD 141569 Imaged with 
NICMOS. Astrophys. J. 525, L53-L56. 

Weinberger, A. J., R. M. Rich, E. E. Becklin, B. Zuckerman, and K. Matthews 2000. Stehar 
Companions and the Age of HD 141569 and Its Circumstellar Disk. Astrophys. J. 544, 
937-943. 

Wetherill, G. W. 1980. Formation of the terrestrial planets. Annu. Rev. Astron. Astrophys. 18, 
77-113. 

Whitmire, D. P., J. J. Matese, L. Criswell, and S. Mikkola 1998. Habitable Planet Formation in 
Binary Star Systems. Icarus^ 132, 196-203. 

Yabushita, S., I. Hasegawa, and K. Kobayashi 1982. The stellar perturbations of orbital elements 
of long-period comets. Mon. Not. R. Astron. Soc. 200, 661-671. 



This preprint was prepared with the AAS lATgX macros v5.0. 



24 



Fig. 1. — Encounter configuration in the frame centered at the primary star with mass Mi. The 
orbit of the passing star with mass M2 is characterized by pericenter distance D, eccentricity e^, 
inchnation i* and argument of perihehon u;*. If length and mass are scaled by D and Mi, the 
encounter parameter are M^, ( = M2/M1), e^, i* and to^. 

Fig. 2. — Left and center panels are the time evolution of orbital eccentricity e and inclination i 
of particles as a function of scaled initial semimajor axis a/D. Time proceeds from (a) top to (c) 
bottom panels. The right panels are face-on snapshots of the disk particles (small dots) and the 
passing star (filled circle). 

Fig. 3. — (a)Orbital eccentricity e and inclination i of particles pumped-up by a passing star, as 
a function of initial scaled semimajor axis a/D, in the case with i* = 5°, w* = 90°, e* = 1 and 
M* = 1. (b)The results with e* = 1 and M* = 0.2 (i* and uj^ are the same). (c)The result with 
e* = 5 and M* = 1. Dashed lines with triangles express analytical expression given by Eqs. (46) 
to (49) in Appendix. The dashed lines (in particular, for i) are almost indistinguishable from the 
numerical results in the inner region. 

Fig. 4. — Orbital eccentricity e and inclination i of particles pumped-up by a passing star, as a 
function of scaled semimajor axis a/D, in the case with lj* = 0°, e* = 1 and M* = 1. Orbital 
inclination of the passing star is (a) 5°, (b) 30°, (c) 45°, (d) 85° and (e) 150°. Dashed lines with 
triangle express analytical expression given by Eqs. (15) and (16) in Section 5. The dashed lines 
are almost distinguishable from the numerical results. 

Fig. 5. — Numerically and analytically calculated cq are plotted, eo is defined by e = (157r/64) 
Co {a/D)^'"^ in the inner disk, in the case with e* = 1 and M^, = 1. (a)The contours of cq 
numerical calculated at every 10 degrees of i^= and a;^,. (b)The contours of analytically calculated 
cq: [cos^a;*{l - (5/4) sin^ i*}^ + sin^w* cos^ i*{l - (15/4) sin^ i^,}^]^/^ (Eq. (15)), as a function of 
i^ and u^. 

Fig. 6. — Numerically and analytically calculated ig are plotted, where io is defined by i = 37r/16 
iQ (a/D)^'"^ in the inner disk, in the case with e* = 1 and M* = 1. Circles are numerical results. 
We calculated the case with oj* at every 10° from 0° to 90°, for each i^,. The results overlap each 
other: iq is almost completely independent of a;*. Dashed lines express analytical expression oi {q: 
sin2i* (see Eq. (16)). 

Fig. 7. — The relative change in semimajor axis Aa after a stellar encounter as a function of initial 
semimajor axis a in the case with i^, = 5°, uj^, = 90°, M* = 1 and e* = 1. 

Fig. 8. — (a) orbital eccentricity e, (b) inclination i, (c) longitude of pericentre uj and (d) longitude 
of ascending node cp of planetesimals, as a function of initial semimajor axis a/D, in the case with 
initial e and i, i* = 5°, lo^ = 0°, e* = 1 and M* = 1. 
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Fig. 9. — (a)ecfr and i^fi, (b)erci and ij-cl and (c)e and i, of the colliding planetesimals, as a function 
of initial seniimajor axis a/D, where we define a is the radius of the colliding point, (cq, io) = (0, 
0) and (i*,u;^,M*,e*) = (5°,0°,1,1). 

Fig. 10. — The same as Fig. 8 except for cq and io are 1 x 10~^. 

Fig. 11. — The same as Fig. 9 except for eg and io are 1 x 10"*^. 
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